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Abstract. Pairs ofnxn matrices whose commutator differ from the identity 
by a matrix of rank r are used to construct bispectral differential operators 
with r X r matrix coefficients satisfying the Lax equations of the Matrix KP 
hierarchy. Moreover, the bispectral involution on these operators has dynami- 
cal significance for the spin Calogero particles system whose phase space such 
pairs represent. In the case r = 1, this reproduces well-known results of Wilson 
and others from the 1990's relating (spinless) Calogero-Moser systems to the 
bispectrality of (scalar) differential operators. This new class of pairs (L, A) 
of bispectral matrix differential operators is different than those previously 
studied in that L acts from the left, but A from the right on a common r X r 
eigenmatrix. 



1. Introduction 

1.1. Background. Let 

(1) CM„ = {(A, Z) £ M nxn | rank ([A, Z] - I) = 1} 

be the set of pairs of complex n x n matrices whose commutator differs from the 
identity by a matrix of rank one. This space arises naturally in the study of the 
integrable Calogero-Moser-Sutherland particle system [vDVOO, Pol06]. In particu- 
lar, the eigenvalues of the time dependent matrix X + itZ 1 ^ 1 move according to 
the ith Hamiltonian of this integrable hierarchy and even allows the continuation 
of the dynamics through collisions [KKS78, Wil98]. 

The KP hierarchy is the collection of nonlinear partial differential equations 

(2) i = 1,2,3.... 

for a monic pseudo-differential operator C of order one whose coefficients are scalar 
functions depending on the time variables ti [SS83, SW85]. If the coefficients 
of C arc further assumed to be rational functions of t\ which vanish as i\ — > 
oo, then the solutions can be written in terms of the matrices in CM„ and the 
poles move according to the dynamics of the Calogero-Moser-Sutherland system 
[AMM77, Kri79, Shi94, Wil98]. This was interpreted as a special case of a more 
general relationship between "rank one conditions" and the KP hierarchy in [GK06] . 

Although it seems at first to be quite different in nature, having no obvious 
dynamical interpretation, the bispectral problem [HK98] turns out to be another 
aspect of this relationship between the KP hierarchy and the Calogero-Moser- 
Sutherland particle system. As originally formulated in [DG86], the bispectral 
problem seeks to find scalar coefficient ordinary differential operators L and A in 
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the variables x and z respectively such that there is a common eigenfunction ip(x, z) 
satisfying the eigenvalue equations 

(3) Lt/j=p(z)4> AV> = 7r(a:)V> 

for non-constant functions p and ir. As it turns out, if one additionally requires 
the operator L to commute with another ordinary differential operator of relatively 
prime order, the solutions to the bispectral problem are exactly the same as the 
rational solutions to the KP hierarchy mentioned above [Wil93]. (Specifically, up to 
trivial rcnormalizations, the bispectral operators are the ordinary differential oper- 
ators that commute with the pseudo-differential operator £ with the identification 
x = t\.) Moreover, the bispectral property for these operators is a manifestation of 
the involution on CM„ given by (X, Z) i— > (Z T ,X T ) which linearizes the dynamics 
of the particle system [Kas95, Wil98] . 

In [Wil98] and the conference proceedings [WilOO], Wilson suggests that the 
correspondence should generalize naturally to the case in which the "rank one con- 
dition" (1) is replaced with a "rank r condition". Indeed, various authors have 
demonstrated that a similar relationship exists between the matrices whose com- 
mutator differs from the identity by a matrix of rank r (1 < r < n), the "spin 
generalization" of the integrablc particle system, and matrix generalizations of the 
KP hierarchy. In particular, the spin generalized system [GH84] was shown to be 
related to the matrix KP equation in [KBBT95] and to the multi-component KP 
hierarchy in [BZN07], and rather general rank r conditions were shown to produce 
solutions to the matrix potential KP hierarchy in [DMH07]. None of these, how- 
ever, has specifically addressed the question of whether and how the results relating 
to bispectrality generalize to the matrix case. 

1.2. Outline. Section 1.3 will introduce a version of the bispectral problem in 
which matrix differential operators act on a common eigenmatrix from opposite 
sides. The main result of this paper will be to demonstrate that this formulation 
allows for the generalization of the results on bispectrality to the spin version of 
the particle system and the matrix KP hierarchy. 

Section 2 introduces the generalization of (1) to the case of arbitrary rank and 
relates it to the dynamics of the spin Calogero particle system. Special attention is 
paid to the block decompositions of the associated operators corresponding to the 
generalized eigenspaces of the matrix Z. 

A wave function and pseudo-differential operator are constructed from a choice 
ofnxn matrices satisfying the rank r condition in Section 3. This r x r matrix 
pseudo-differential operator is shown to satisfy the Lax equation of the KP hierarchy 
(2). A key component of the proof is the explicit construction of an rn-dimensional 
space of finitely supported distributions in the spectral parameter which annihilate 
the wave function. 

Several obvious group actions on the space of matrices satisfying the rank r 
condition are investigated in Section 4 with emphasis on their effect on the corre- 
sponding KP solution. Of special interest is the bispectral involution which has the 
effect of exchanging variables and transposing the wave function. 

The main theorem is the construction in Section 5 of commutative rings of matrix 
differential operators in x and z and the demonstration that they have the KP wave 
function as a common eigenfunction. 
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A final section contains closing remarks and lists problems for future research 
on this topic. 

1.3. Notation and Matrix Bispectrality. We will make use of the notation M L 
and M R to distinguish between the cases in which the operator M is acting from the 
left or the right, respectively. So, for instance, if M and P are both r x r matrices, 
then 

[M,P] = (M L -M B )P. 

Similarly, if 

N 

(4) L = J2M l (x)dl 

i=0 

is an ordinary differential operator in x of degree N with coefficients Mi(x) that 
are r x r matrices and ip(x) is an r x r matrix function we define 

MVO = Lty) = J2 Mi(x) ( JU(*)) , 
as usual. However, the operator can also act from the right 

N / fji \ 

w) = E \w^ x) ) Mi{x) - 

i=0 ^ ' 

Equivalently, if we denote by L T the differential operator with coefficients Mj that 
are the ordinary matrix transpose of the coefficients of L, we can say 

(5) L R (V0 = (£ T W< T )) T . 

Definition 1.1. A bispectral triple (L,A,if>) consists of a differential operator L 
in x as in (4), a differential operator A in the variable z also having r x r matrix 
coefficients, and anrxr matrix function ip(x, z) of x and z satisfying the equations 

(6) L L (tp) = p(z)ip and A K (ip) = ir(x)ip, 
where p(z) and ir(x) are non-constant, scalar eigenvalues. 

This seems to be a natural matrix generalization of the scalar bispectral problem 
for differential operators considered in [DG86]. However, we note that this differs 
from the matrix generalization previously considered by Zubclli [Zub90] in which 
both operators acted from the same side, and also from the "bundle bispectrality" 
considered by Sakhnovich-Zubclli [SZ01] where the operators L and A were allowed 
to depend on both variables. (Here we are interested only in the case that L is 
independent of z and A is independent of x.) 

In the rest of the paper we will use the notation Ik for the k x k identitity 
matrix. Also we will abuse notation by using the symbol I to denote the identity 
transformation on many different vector spaces whenever its use should make it 
clear which is intended. 
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2. Spin Calogero Matrices 

Let sCM^ be the the set of 4-tuples of matrices (X, Z, A, B) such that the n x n 
commutator [X, Z] differs from the identity by the rank r matrix BA: 

(7) sCM r n = {(X,Z,A,B) \X,ZeM nxn , A,B T E M rxn , [X, Z]—I = BA ^ 0}. 

This space arises naturally in the description of generic initial conditions of the Spin 
Calogero particles, as we will see below. More importantly, the dynamics linearizes 
there (when one considers the phase space to be sCM^ modulo the action of GL(n) 
to be described in the section on symmetries). 

Let qi (1 < i < n) be the distinct positions of n particles on the complex plane, 
qi be their momenta, and = (3jCti be their "spins" represented as the products 
of n column r- vectors and n row r- vectors (3j subject to the constraint fa = — 1. 

We associate to this data the matrices (X, Z, A, B) e sCM^ in the form 



Xij — qiSij Zij — qiSij + (1 — Sij) 



fij 



Qi - <lj 

(Pi 

A = (ai • • • a n ) B = : 

The dynamics of the eigenvalues of X + itZ l ~ 1 are governed by the Hamiltonian 
H = trZ 1 . This is the spin Calogero system [GH84, Kri79]. In the special case r = 
1, this reduces to the more famous (spinless) Calogero-Moser-Sutherland particle 
system [vDVOO]. 

2.1. Block Decomposition. For a fixed choice of (X, Z, A, B) G sCM^ we get a 
decomposition of V = C™ into generalized eigenspaces of Z: 

V = v * V\ = {v &C n \(Z — XI) k v = for some k > 0}. 

A 

The restriction of Z to V\ will be denoted by 

Z x = XI + N x 

where / is the identity operator on V\ and N\ is nilpotent. Rational expressions 
in zl — Z\ below will always be interpreted by expansion in positive powers of the 
nilpotent N\. 

We will utilize subscripts A and \i which will run over the eigenvalues of Z to 
similarly denote the blocks of other linear operators associated to this decomposition 
of C™. Specifically, A\ : V\ — ► C will be the restriction of the map A, B\ : C r -» V\ 
will be the map B followed by projection onto V\, and for a linear operator M from 
C to itself (such as X) will be the block corresponding to the map from 

to V\. 

The sCM condition (7) involves the commutator [X, Z] = (Z R — Z L )(X). Interest- 
ingly, although the operator Z R — Z L = —&d(Z) is not invertible, its "off-diagonal" 
action is invertible which allows us to solve for X\^ when \i ^ A. 

Lemma 2.1: Let A ^ /x be generalized eigenvalues of Z. Then 

X* = ((Z A ) E - (Z„) J" 1 B,A X = g ^^ BA 
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Proof: Since 

(Z K - Z L )^ = (A - m)/^ + (N x \ - (A^) L 
differs from the a nonzero multiple of the identity by a nilpotent matrix, it is 
invertible. Specifically we may invert it in general using 



(7 WJ* 



fc>0 

Applying this when solving [X, Z] — I = BA for any off-diagonal block of X yields 
the claimed formula. □ 

It will later be necessary to evaluate residues of matrix functions written in terms 
of the blocks Z\. For this purpose the following "obvious" lemma will be useful. 
For convenience we introduce notation for "divided derivatives" : 

k\ dz k 

Lemma 2.2: Let f(z) be a rational function that is regular at z — A, then 
Proof: 



z =x \ (zl - z x ) k+1 J z=\ \ J w M zI-Zx 



7=x V y ~' zl - Z X J "=>• \ J K ~' ^ (z - A) s+1 



= Res I fW( z )J2(-d z y/ S \-^- I = Res 

z=X \ ^ (Z — \) ] z=X 



Y,d s J S \.f [k] {z)N{ 



s>0 



- A 



Res (fl% + N x )-!—) =fW(Z> 
z=x y z — a ) 



)■ 



□ 

Of course, in the above lemma f^ k \Zx) is a matrix and may not commute with 
other matrices appearing. So, one needs a little care in applying the Lemma 2.2. 

3. Matrix KP Hierarchy 

Let v = (X, Z, A, B) G sCM^ and associate to it the wave function tp v depending 
on the spectral parameter z and the times t = (t\, £2, h, • • -) : 

(8) Mt, z) = 7(*i *) + AX-\zI n - Z^B) , 

where 1 



00 

,-- 

i=l \i=l 



X = = iUZ 1 - 1 - X, and 7$ z) = exp ( 



^^Note that the dependence on t{ in —X is such that its eigenvalue dynamics are governed by 
the i th spin Calogcro Hamiltonian. 
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If q v (z) = det(zI-Z) is the characteristic polynomial of Z then the wave function 

(8) can be multiplied by q v and an exponential so as to yield a polynomial in z 
with coefficients that are rational in the times 

(9) K(t,z)^ 1 - 1 (t,z)^(t,z)q,(z). 

Indeed, q v (z)(zl — Z)^ 1 is the classical adjoint of zl n — Z is and hence polynomial in 
z. Letting d = be the differential operator in x = ti, we note that the ordinary 
differential operator K v = K(t, d) satisfies 

(10) ^ v (t iZ ) = J—r^iz). 

qv\z) 

The main goal of this section is to prove that the pseudo-differential operator 

L v = K v odo jC^ 1 

is a solution to the matrix KP hierarchy in that it satisfies the Lax equation (2). 
As in [Wil93] (see also [Kas95, SW85]), the proof will involve identifying finitely 
supported distributions in z that annihilate the function 

3.1. Conditions satisfied by ip v . Consider a generalized eigenvalue A of Z with 
multiplicity I and use the notation of Section 2.1 to denote by Z\, Xx^, Ax, etc. 
the blocks of the operators X, Z, A and B. Let v 6 C r£+£ have the decomposition 

/ v ° \ 

Vl 

v 2 

V = 

ve-i 
\w J 

where Vi £ C r and w <E C £ and define the distribution c Vt x taking r x r matrix 
functions of z to r component constant vectors by the formula 

(11) c v , x (f(z)) - Res (f(z) ■ [Ax(zl - Z^w + ^> - \) l vA ) . 



i=0 



In this section we will show that there are rl linearly independent distributions of 
this form satisfying c v ^x(ipv{t, z)) = 0. Consequently, by running through all of the 
eigenvalues of Z we obtain in this manner an rn-dimensional space of conditions 
satisfied by the wave function. Indeed, if {Ai} are the generalized eigenvalues of 
the n x n matrix Z with multiplicities {li}, then n — Y^U- 
Consider the r x (rl + £) matrix 



Tx = (B x N x Bx NlBx 



Nt'Bx 



-Xxx) ■ 



Lemma 3.1: If v G kcrTA then c v ^x (Vv(£, z)) = for all values of the variables t. 
Proof: Note first that ipv{t, z) has the block decomposition 



(12) 



Mt, z) = j(z, t)(i-J2 A^x- 1 )^ - Z^BA 



where again the sum is taken over all (not necessarily distinct) pairs of generalized 
eigenvalues k and \i of Z. 
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Now, we wish to use Lemma 2.2 to expand the residue in (11) where f(z) is 
replaced by (12). It will be convenient to introduce the abbreviation 

Cfi — — ^ A K X K X , 

K 

so that we have 

(13) z) = 7$ z ) ( 7 + E c ^ zI - z ^y lB ^j 

and 

(14) A x = -Y,CpX,>x. 

The various contributions to the residue are usefully organized according to 
dependence on C^. First of all, there is the contribution independent of C^. It is 
given by 

(A) Res (7$ z)A x (zI - Z x )~ l w) = A xl {t, Z x )w. 

Making use of Lemma 2.1 one finds that the contributions containing for \i ^ A 
are 



(B) Jff i( zI - Z^r l B^A x {zI - Z^w^t, z)) = C^X^t, Z x )w. 
Next we turn to the terms involving C x . The first one is 

C\ Res ^f(t, z)(zl - Zx^Bx J2(z - \yv^j 

= c„ h& z )( zI - z x y +1 N{B x (z - xy Vj ) 

i,j=o 2 

t-i 

(C) =C xl (t,Z x )J2N l x B x v i . 

i=0 

The other term linear in C x is 

C\ Res ( 7 (t, z)(zl - Z x )- l B x A x {zI - Z x )- X w) 
=C X Res (76 z)(zl - Z X )-\[X XX , Z x ] - I)(zl - Z x )~ l W ) 
= -C x Res (<y(t, z)(zl - Z x y 2 w) 

+ C\ Res (7$ z)(zl - Z X )-\[X XX , Z x ] - z)(zl - Z^w) 
= - C xl '(t, Z x )w - C x Rcs( 7 (C z)(zl - Z^X^w) 

+ C x Res(X X x(zI - ZJ-^it, z)w) 

z=X 

= - Cxi {I Zx)w - C xl (t, Z x )X xx w + C x X xxl (t, Z x )w 

(D) =C x X xx7 (t, Z x )w - C xl (t, Z x )X xx w. 
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Since v 6 ker T\ is equivalent to the statement 

(15) J2 N * B ^- X ^ w = > 

we see that (C) cancels against the second term in (D). So, combining all four 
terms gives 

(A) + (B) + (C) + (D) = A xl (t, Z x )w + ]T C^X^it, Z x )w + C x X xxl {t, Z x )w. 

Applying (14) shows that this is equal to zero as required. □ 

Lemma 3.2: The distributions c VtX for v 6 kcrl\ form an ^-dimensional space. 

Proof: Note that the map il : v 6 kerl\ i— > c V:X is itself a linear map. What we 
need to prove, therefore is that 

dim kcr T x — dim ker £1 = r£. 

A vector v clearly does not lie in the kernel of £1 if Vi ^0 for any i. The dimension 
of the kernel of il is therefore equal to the dimension of the space of vectors w with 
the property that X xx w = and A x (zl — Z x )^ 1 w = 0. In fact, we will show that 
the only such w is the zero vector (and hence that dim ker il = 0). Beginning with 
the fact that 

[(zI-Z x ),X xx ]-B x A x =I. 

Multiplying by (zl — Z x )~ x on the right, applying both sides of the resulting equa- 
tion to w and then multiplying by (zl — Z x )~ x on the left gives us that 

X xx (zl - Z x )- l w = {zl - Z x y 2 w. 

Expanding both sides of this equation in terms of powers of (z — A) and equating 
like powers gives us that 

X xx N^w = kN^w, for k> 0. 

Since A^a is nilpotent, for a sufficiently large k the left-hand side is equal to zero. 
But the equation then tells us that N^~ lr w is then also equal to zero, which again 
means that the left hand side would be zero for a smaller value of k. Repeating 
this process until k = 1 we find that w — 0. 

A similar argument shows that dim ker T x = r£. Considering instead the vectors 
w such that w t Fa = implies that w T X = w 1 (zl — Z X )~ X B X = and the same 
process reveals that w — so that T x has rank t. Consequently, its kernel has 
dimension (ri + £)—£ = r£. □ 

3.2. The Kernel of Id,. The results of the previous section on the distributions 
annihilating ip v give us information about the kernel of the matrix ordinary differ- 
ential operator /C„ defined in (10): 

Corollary 3.3: Let c VjX be as in Lemma 3.1. Then the r component vector valued 
function 

(16) 

is in the kernel of the operator K v . 
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Proof: Since c Vt \ commutes with multiplication and differentiation in x = ti, we 
have 

)C v (f> v ,\ = c„ ; a (K. v ^(t,z)q~ 1 (z)) = c Vt \(^ v ) = 0, 
by (8) and Lemma 3.f . □ 

In fact, the entire kernel of /C„ is spanned by functions of this form, and as a 
consequence they satisfy certain useful linear differential equations. 

Corollary 3.4: If (j>(i) is a vector in the kernel of Id, then it is a linear combination 
of the 4> v ,\{i), and so satisfies the equation 

Qk Q 

^ = w k *®- 

Proof: By making use of all of the eigenvalues of Z, Corollary 3.3 gives us rn 
linearly independent vector functions in the kernel of the nth order ordinary differ- 
ential operator with r x r matrix coefficients. Since they are linearly independent 
(those corresponding to the same eigenvalue are linearly independent by Lemma 3.2 
and those corresponding to different eigenvalues cannot be linearly dependent due 
to the factor of e xX ) this accounts for the entire kernel of K v . 
Note that j(t, z) trivially satisfies 

Qk Q 

^M=AM = ^7(t;*). 

Now, the proof here is elementary because differentiation in ti commutes with the 
residue, multiplication by functions of z and matrix multiplication in the definition 
of <fiv,\ and applies by linearity to the entire kernel. □ 

3.3. The Lax Equation. Now we come to main point of this section. If the v 
moves according the spin Calogero dynamics the wave function ip v depends on the 
time variables t, and this produces a solution of the matrix KP hierarchy. More 
precisely: 

Theorem 3.5: The pseudo-differential operator C v = K v o d o K,~ x satisfies 

d 

gf^ v = hi- 
proof: First, we note that the (pseudo)-differential operator {L % v )-oK, v is actually 
a differential operator since 

(4)_oAC„ + {C v )+oK v = K v o d* 

and therefore o JC V = — {C l v )+ ° K, v + K, v o d l . 

Now, let <p(x) be a vector function in the kernel of the operator JC U . Then 
applying J^- to the equality /C„0 = and using Corollary 3.4 we find 

= ^o£„(0) = 

= {K v ) u( j> + C v o }C V (4>) 

= (£„) t4 (0) + (4)+ o K v {<j>) + o JCM) 
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However, since <fr is in the kernel of K, v we know that (C\,)+ ° K. v {4>) = 0. Then the 
last displayed equality gives us that the entire kernel of K, u is in the kernel of the 
ordinary differential operator (K, v ) ti + {£-{,)- fcv Since this operator has order 
strictly less than n, it can only have such a large kernel if it is the zero operator 
and we conclude 

{K v ) ti = -(£*)_ o/C„. 

Using this we find that 

(c v ) tt = (q^so^-^oflo^-'o^or 1 

= -(Cl)- o K v o d o k- 1 + K v o a o /C; 1 o (4)- o K v o k- 1 

= [£„,(4)_] = [(4) + , J c w ]. 

□ 



4. Symmetries 

The symmetry Ih1| cZ j of sCM^ induces the integrable dynamics of both 
the particle system and the wave equations of the matrix KP hierarchy. Here are 
some other symmetries and how they affect the KP Lax operator. 

4.1. Action of GL(n). For any G G GL(n) we define Sg acting on sCM^ by 

Sg : sCM; -» sCM; 

(X,Z,A,B) i ^ (GXG- 1 ,GZG- x i AG- 1 ,GB). 

The Lax operator £ v is unaffected by this action: Cs G {v) — ^-u ■ 

Since this symmetry docs not affect the corresponding dynamical objects from 
the previous sections, it makes sense to consider sCM^ modulo this group action as 
the phase space of the spin Calogero particle dynamics as well as the corresponding 
matrix KP solutions. 

4.2. Action of GL(r). Clearly, if we have G G GL(r) then we can conjugate 
solutions to the matrix KP hierarchy to get solutions that are technically different, 
but not very different. This also manifests itself as a group action on the level of 
sCM^. Let G G GL(r) then if 

s G : sCM; -> sCM; 
(X,Z,A,B) i ► (X,Z,GA,BG- X ) 

one finds that C- Sa {v) — GC V G~ X . 

4.3. Changing r. There is an easy way to take an r x r solution and turn it into 
an R x R solution for r < R. Let a be an R x r and b an r x R matrix such that 
ba = I is the r x r identity matrix. Defining U a .b : sCM^j — ► sCM^ by 

U aib (X, Z,A,B) = (X, Z,aA,Bb). 

one then has 

£u a , b (v){t) = aC v (j)b. 



SPIN CALOGERO AND BISPECTRAL MATRIX KP 



11 



4.4. The Bispectral Involution. Finally, we have an important discrete symme- 
try whose effect on the KP solution will be the subject of the next section: 

b:sCM; -> sCM r n 
v={X,Z,A,B) i ► z> = (Z T ,X T ,B T ,A T ) 

The significance of this symmetry on the KP solution is most easily seen by 
looking at the wave function ip v as a function of x = ti and z only (setting all of 
the other times equal to zero). As in the case r — 1, it involves an exchange of x 
and z, but when r > 1 one must also take the transpose of the function: 

(17) i/> v (x,z)=il%(z,x). 

5. BlSPECTRALITY 

Let v G sCM^ and q v (z) — det(zl — Z). In this section, since the dynamics are 
not significant, we will consider x = t\ and ti = for i > 1. Thus, for instance, we 
will write tj} v {x, z) for ip v ((x, 0, 0, . . .), z) and 7(2;, z) = e xz . In the next section we 
will associate two commutative rings of ordinary differential operators (one acting 
in x and one acting in z) to the choice of v G sCM^ and then in the following section 
we will demonstrate that ip„(x, z) is a common eigenfunction for the operators in 
the rings. In particular, any operator from each of the rings along with ip v form a 
bispectral triple as in Definition 1.1. 

5.1. Commutative Rings of Matrix Differential Operators. Associate to 
v G sCM^ the ring R v C C[z] defined by the property that the polynomials preserve 
the conditions annihilating ip v (t,z) from Lemma 3.1: 

R v = {p G C[z] I c„ ;A {tpu) = => c ViX (p(z)il> v ) = 0} . 

Lemma 5.1: The ring i?„ is non-empty. In particular, 

ql(z)C[z]cR v . 

Proof: Note that q v (z)ip u (t, z) — K, u ~f(t,z) is non-singular in z. Then the claim 
follows from the fact that c v ^\(q v (z)f{z)) = for any non-singular function /. □ 

By substituting the pseudo-differential operator L v into these polynomials, we 
associate a commutative ring of pseudo-differential operators 

Tl v = {p{C v )\p£R v } 

to v. However, as the next lemma demonstrates, these are in fact differential 
operators. 

Lemma 5.2: If p G R v then L = p{C v ) is a differential operator (as opposed to a 
general pseudo-differential operator) satisfying the eigenvalue equation 

Lip v (x,z) =p(z)ip v (x,z). 

Proof: Since the leading coefficient of /C„ is a nonsingular matrix (in fact, it 
is the identity matrix because of the form of ip v ), it is sufficient to show that the 
kernel of /C„ is contained in the kernel of /C„ o p(d) because then we know that 
this ordinary differential operator factors as L o /C„ for some ordinary differential 
operator L which meets all of the other criteria. 
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So, now let <j)(x) be a function in the kernel of /C„. By Lemma 3.3 we know that 
4>(x) is a linear combination of functions of the form (16). However, 

K v op{d)c v ^ x {l{x,z)q~ 1 {z)) = c VtX {p{z)ip u {x,z)) = 

by (10) and the definition of R v . □ 

We will also associate a commutative ring of ordinary differential operators in z 
to v. Applying the procedure above to the point v° = (Z T , X T , B T , A T ) E sCM^ 
we have another commutative ring 1Z V \> of differential operators in x. We convert 
them to differential operators in z by simply replacing x with z, d x with d z and 
transposing the coefficients: 

Ki = {L T (z,d z )\L(x,d x ) £1Z V \>}. 

5.2. Common Eigenfunction. Our main result is the observation that ip v (x,z) 
is a common eigenfunction for the differential operators in the rings 1Z U and Vj v 
satisfying eigenvalue equations of the form (6) : 

Theorem 5.3: Let p E R v and ir E R v \>, then there exist ordinary differential 
operators L(x,d x ) € 1Z U and A(z,d z ) G IZ^ such that 

Lip„(x,z) =p(z)ip v (x,z) and A R ip„(x, z) = ir(x)ip„(x, z). 

Proof: The first equation follows from Lemma 5.2. Similarly, it follows from 
Lemma 5.2 that there is a differential operator Q(x,d x ) with the property that 

QVV {x, z) = n(z)4>„t, (x, z). 

Exchanging the roles of x and z in this equation, taking the transpose (see ), and 
applying (5) and (17) results in the second equation of the claim. 

□ 



6. Example 

For the sake of clarity, we briefly illustrate the main ideas with an example. 
Consider v = {X, Z, A, B) E sCMj where 



X = -1 -1 z = 




1 
1 



B 





1 





(• 



















1 N 








) 


I— 1 







Then 



/ / -2zx' 1 +3zx+2x-2 1 \ \ 

Mx,z) = e"[l+[ {x - x 2 z ] -2 z2 ~~' ^W z ) ■ 

\ \ (x-2)xz 2 (x-2)xz J ) 

This can be written as ip v = K v e xz / ! q u {z) where q v {z) = det(z/ — Z) = z 3 and 



IC=[^W Jd+ W K+ n i P 3 

\ _ (x-2)x U / V x-2 (x-2)x ) V U V 

is an ordinary differential operator. 
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-1 














-1 






To find the conditions satisfied by ip v {x, z), we note that the kernel of 

r = 

is made up of vectors of the form 

v = (ci c 2 c 4 c 5 c 6 3ci+c 2 c 8 -ci-c 2 ) T . 

Hence, we conclude that c v fi{^ v {x, z)) — where 

««(/(»))- a 

Moreover, c v fi(p(z)ij} u ) — whenever p £ z 4 C[z] — Consequently, we can 
find an ordinary differential operator having any of these polynomials as its eigen- 
value. In particular, solving 

K v o d 4 = L o K, v 

for L wc find 

( 8(3x 4 -18x 3 +50x 2 -63x+30) 4(23x 3 -93x 2 + 138x-72) 

(x-2) 4 x 4 (x-2) 4 x 5 
8(2x 3 -5x 2 +9x-6) 4 (2x 4 -8x 3 +27x 2 -42x+24) 

(x-2) 4 x 3 (x-2) 4 x 4 

(4(6x 3 -27x 2 +49x-30) 4(l3x 2 -35x+26) \ 

(x-2) 3 x 3 (x-2) 3 x 4 j q 

4(4x 2 -7x+6) 4 (2x 3 -6x 2 + 13x-10) I 
(x-2) 3 x a (x-2) 3 x 3 / 

( 8(x 2 -3x+3) 4(3 x-4) \ 

(*-2) 2 x- M 5 ? x «2 , d 4 

4 4(x 2 -2x+2) I + ° 
(x-2) a (x-2) 2 x a / 

which satisfies Li/v = z 4 Vv 

Of course, we can follow this same procedure beginning with another element of 
sCMg. In particular, if we begin with 

is b = {Z T ,I T ,B T \A T ) 

instead then the differential operator we produce will be 

4(l6x 2 +65x+90) go x+ 21§ \ / 4(l6x 2 +65x+9p) 8(10 x+27) 

8x + 12 4(4x 2 +21x+36) + I 4(2x+3) 4(4x 2 +21 x+36) 



Q 



7T 



/ 2(l2x 2 +65x+90) 40 x+108 \ / 34x+60 4(2x+9) \ 

/a 12 6 \ 

+ ^ Q ^ 4 _f 6 Ja 4 -49 5 + 9 6 . 

The function if)^ is an eigenfunction for this operator satisfying QVv i x , z ) — (4z 4 — 
4z 5 + z 6 )Vvb(a;,z). 

More interestingly, since ip v *(x, z) — ipj(z,x), if we transpose the matrix coeffi- 
cients on Q, replace x with z and d — d x with d z , we get a differential operator A 
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in the variable z. This operator applied to ip v (x,z) (the wave function computed 
earlier) from the right satisfies 

A R ip u (x, z) = (Ax 4 - Ax 5 + x & )ip u (x, z), 

demonstrating bispectrality. 

7. Conclusions and Comments 

The main results of the present paper can be viewed as another step in address- 
ing the "bispectral problem" of F.A. Griinbaum [HK98, DG86], seeking operators 
satisfying eigenvalue equations of the form (3). In [DG05], the authors considered 
the case in which one of the operators is a second order difference operator with 
matrix coefficients and the two operators act on matrix eigenfunctions from dif- 
ferent directions. However, bispectrality for matrix differential operators has only 
been studied with both operators acting from the left [Zub90]. Here we consider 
the case (6) in which the operators are r x r matrix differential operators acting 
from different directions. Since our construction conveniently reproduces the re- 
sults of Wilson's seminal paper [Wil93] in the special case r = 1, this particular 
formulation of the bispectral problem appears to be the correct one for generalizing 
those results to the case of matrix differential operators. However, the method of 
proof and especially the explicit formulation of the "conditions" satisfied by the 
wave function above are novel even for r = 1. 

In [Wil93] , it was shown that the bispectral operators associated to sCM* are in 
fact the only bispectral scalar ordinary differential operators which commute with 
operators of relatively prime order up to obvious renormalizations and changes of 
variable 2 . By Lemma 5.1 it follows that the differential operators produced by the 
construction in this paper also all have the property that they commute with other 
differential operators of relatively prime order. 

In addition, this paper can be seen as contributing to the literature establishing a 
link between bispectrality and duality in classical and quantum integrable systems. 
(See, for instance, [FGNROO, Hai07, Kas95, KasOO, KasOI, Wil98].) Again, the 
main results of the present paper for the spin Calogero system in the case r = 1 
reproduce results previously presented for the spinless case in [Kas95, Wil98]. 

Some questions arise naturally which we have not pursued. There arc additional 
commuting Hamiltonians for the spin Calogero system [GH84] and corresponding 
isospectral deformations for the multi-component KP hierarchy [BtK95], but their 
relationship to bispectrality has not been explored here. We have not looked at 
the algebro-geometric implications of the rings TZ V . Certainly as in the case r = I 
[Wil93, Wil93], these contain operators of relatively prime order are isomorphic to 
the coordinate rings of rational curves with only cuspidal singularities. However, 
whether there is any further algebro-geometric significance such as was found in 
[BW02] or whether every commutative ring of matrix ordinary differential opera- 
tors with these properties is necessarily bispectral have not been considered. These 
questions, along with the obvious question of what other matrix differential opera- 
tors satisfy equations of the form (6) will hopefully be addressed in future papers. 
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These were called "rank one" operators in that context, but we will avoid that terminology 
here to avoid confusion with the rank r which is something different. 
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